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Abstract 

We consider the anomalous sub-diffusion of a class of Gaussian pro- 
cesses that can be expressed in terms of sums of Ornstein-Uhlenbeck 
processes. As a generic class of processes, we introduce a single param- 
eter such that for any v G (0, 1) the process can be tuned to produce 
a mean-squared displacement with E[a; 2 (t)] ~ t v for large t. 

The motivation for the specific structure of these sums of OU pro- 
cesses comes from the Rouse chain model from polymer kinetic theory. 
We generalize the model by studying the general dynamics of individ- 
ual particles in networks of thermally fluctuating beads connected by 
Hookean springs. Such a set-up is similar to the study of Kac-Zwanzig 
heat bath models. Whereas the existing heat bath literature places its 
assumptions on the spectrum of the Laplacian matrix associated to the 
spring connection graph, we study explicit graph structures. In this 
setting we prove a notion of universality for the Rouse chain's well- 
known E[a; 2 (i)] ~ scaling behavior. Subsequently we demonstrate 
the existence of other anomalous behavior by changing the dimension of 
the connection graph or by allowing repulsive forces among the beads. 

1 Introduction 



Due to recent and compelling experimental observations using advanced mi- 
croscopy [H El [3] there is theoretical interest [H H3 [HJ [91 [10] in anomalous 
diffusion - stochastic processes whose long-term mean-squared displacement 
(MSD) satisfies E[x 2 (t)] ~ t v where v ^ 1. In each of the cited references, 
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the observed behavior is sub-diffusive, where v £ (0, 1). The canonical ex- 
ample of a sub-diffusive process is fractional Brownian motion [8], but in 
this paper, we focus on a model from polymer kinetic theory and natural 
generalizations. 

The Rouse chain model of a polymer is a series of thermally fluctuat- 
ing beads {x n (t)}t>o, n G {1, . . . , N} that interact with nearest neighbors 
through linear spring forces. It is a standard observation in the physics 
literature [TTJ [121 E] that while the center-of-mass of the chain is a diffu- 
sive process, there exist positive times t\ and tn such that individual beads 
roughly exhibit the following MSD profile: 



E[*n(*)] 



t, t < n 

n<.t<.T N (l.i) 

t, t » Tjv 



The times ri and tat are called the first and last relaxation times respectively. 
Two of the primary projects in the present paper are to give precise mathe- 
matical meaning to a profile such as (jl.ip and to show that the anomalous 
exponent v = ^ in the intermediate timescale is determined by the geometric 
structure of the graph of connections among the beads. 

The observation that the dynamics of an individual particle in a network, 
a so-called distinguished particle process, can exhibit anomalous diffusion is 
not restricted to polymer kinetics. In a series of papers [HI EJ [J_5] the 
authors studied the behavior of a distinguished particle in a Kac-Zwanzig 
heat bath, a model used in molecular dynamics theory to study the force 
exerted on a particle by a randomly fluctuating environment. Of present 
interest are the articles [5] and [8] wherein the authors showed that in an 
appropriately constructed large- N limit, a family of distinguished particle 
processes can converge weakly to a sub-diffusive limiting process. In [8], 
this process is fractional Brownian motion, while in [5] the limiting process 
is the so-called generalized Langevin equation (see also |16| ) with a power 
law memory kernel. It is worth noting that in each of the above cases, the 
results followed from assumptions that were placed on the spectrum of the 
weighted spring connection graph, rather than directly on its weights and 
geometric structure, which is the goal of the present paper. 

In [10], the authors introduce a common mathematical framework to 
address the sub-diffusion seen in these models: a class of Gaussian processes 
expressible in terms of a Brownian motion plus a sum of Ornstein-Uhlenbeck 
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processes, 

N-l 

x(t) = c B {t) + c kz k {t) (1.2) 

k=l 

where each Zk satisfies the SDE 

dz k (t) = -XkZkit) + dB k {t) (1.3) 

The collection of standard Brownian motions {i?fc,A r }fc<Af-i & r e assumed 
to be independent. The positive constants {^k}k<N-i will be called the 
diffusive spectrum and the constants {ck}k<N-i will be called the coefficient 
family. The inverses of the elements of the diffusive spectrum := A^ 1 
are called the relaxation times of the process. Henceforth we will refer to 
processes defined by (jl.2p and (jl.3p as SOU processes. 

One can think of an SOU process as a diffusing particle that is free 
to explore all of space, but is subject to a sequence of linear mechanical 
responses from its environment. The central notion of this paper is that 
anomalous diffusion can arise from the structure of timing of this cascade of 
responses. One can directly show (see Section 1 1 . 2f) that when the beads in 
a network interact through linear spring forces, the associated distinguished 
particle process is exactly expressible in terms of equations (II. 2\\ and fjl .31) . 

In [10], the authors laid out the general relationship between the diffusive 
spectrum {^k}k=i anc l the intermediate timescale anomalous exponent for 
SOU processes. Generalizing the Rouse spectrum A& = sin 2 (/c7r/2A r ) from 
the polymer kinetic theory [HJ [13] by defining the diffusive spectrum to be 

A^=(£)V> (1.4) 

we find that SOU processes can exhibit any desired anomalous exponent 
between and 1. Indeed, the full MSD profile of an SOU process with 
generalized Rouse spectrum is given by 



E[x 2 (t)] 



t, t<^Tl 
t l ~~P, T G <^t<^T N (1.5) 
t, t » Tn 



for any p > 1. The actual exponents observed in experiments vary widely, 
and so the existence of stochastic processes that exhibit robust and var- 
ied anomalous behavior is appealing to experimentalists and engineers. At 
present, this community still lacks effective statistical inference tools for 
these sub-diffusive processes as well as the ability to conduct simulated ex- 
periments that are not computationally prohibitive. 
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1.1 Summary of Results 

In this paper, we seek to give a rigorous interpretation to the claims made 
in the companion work [10] . In Theorem 12,11 we give precise meaning to the 
MSD profile (|1.1|) . The key observation, independently noted in [8], is that 
for these models, the largest relaxation time rjy tends to infinity with the 
number of modes N. The family of processes is tight, with a sequence of 
autocorrelation (ACF) functions that converge uniformly on compact sets. 
Therefore, demonstrating anomalous diffusion for large t in the limiting 
process is tantamount to proving it for the intermediate timescale of the 
finite N processes. 

There is a robust sense in which the anomalous exponent is determined 
by the diffusive spectrum rather than the coefficient family. This fact is vital 
if one hopes to perform statistical inference on an SOU process, see |10j for 
further discussion. In Proposition 12.21 we show that setting the coefficients 
to be i.i.d. random variables with mild restrictions, the exponent v remains 
the same. 

While the generic SOU structure can support all anomalous sub-diffusive 
exponents it turns out that the behavior of distinguished particle processes 
is not so varied. It is conjectured in the physics literature [11] that the 
Rouse model is a universality class that captures the qualitative behavior 
of a wide variety of bead-spring networks. The main result of Section [3] is 
that this is indeed true in some rigorous sense (Theorem 13. 2p . The only 
exponent seen for a wide class of models is v = |. However, it is possible to 
construct weighted networks that produce different behavior by 1) changing 
the dimension of the underlying spring connection graph, Section 13.31 and 
2) by allowing for repulsive forces among the beads, Section [331 Ultimately, 
such modifications cannot account for the wide behavior seen experimentally. 
One will likely need to account for some combination of hydrodynamic self- 
interaction [T7] [18] and excluded volume effects [IT] , but rigorous study 
of these effects without pre-averaging approximations remains an unsolved 
problem. 

1.2 Distinguished particles in bead-spring networks 

We demonstrate the connection between distinguished particle processes and 
SOU processes. Let x(t) = . . . , X2(t), . . . , xjv(i)) denote the 

locations in M. d of a set of particles at time t > 0. For the sake of simplicity 
we take d = 1, although this is not essential (see Section [3.1.2|) . Following 
the development of flexible polymer kinetics [111 I13| . the particles are subject 
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to random thermal fluctuations while interacting through a given quadratic 
configuration potential 

*(x) = -J2Kn m \Xn-Xm\ 2 . 

The set of pairs £ := {(x n ,x m ) : K nm > 0} constitutes the set of edges of 
the graph Q associated with network. 

Particle dynamics are formally set by a balance of forces through the 
Langevin equation, 

mx = -rjx(t) - V*(x(t)) + aW(t) 

where r] is the viscosity of the fluid in which the particles are immersed, and 
the strength of the noise a is related to the viscosity through the fluctuation 
dissipation theorem: a = \J2kBTr]. The constant T is the temperature of 
the fluid and ks is Boltzmann's constant. For simplicity we will renormalize 
the dynamics so that 77 = 1. 

The common mass m of the particles is considered to be small. In the 
companion paper [10], we observe that the zero- mass limit is singular and 
non-trivial to analyze. Here we restrict our attention to the weak (over- 
damped) zero-mass limit, as described in [10], which amounts to setting 
m = and heretofore taking the stochastic forcing term W(t) to be an i.i.d 
vector of standard Brownian motions. 

The force exerted by the configuration potential on the n-th bead is 

^ 'x„^ — ^ ^ l^"nm(%m ■Em) 

leading to the linear system of SDEs 

dx(t) = Lx(t)dt + ad\V(t) (1.6) 

where L is the so-called Laplacian matrix for the spring connection graph 
Q. The Laplacian matrix is sometimes written L = A — D where A is the 
weighted adjacency matrix for Q and D is a diagonal matrix whose entries 
are the sums of spring constants D nn = ^ m n nm . 

We note that L is symmetric and negative definite. As such it can be 
diagonalized in the form 

L = QAQ 1 (1.7) 

where A is a diagonal matrix with the eigenvalues of L as its entries, and Q 
is an orthogonal matrix with the eigenvectors of L as its columns. We make 
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use of a few standard observations. First, is always an eigenvalue of L and 
if Q is connected, then has multiplicity 1. The eigenvector associated to 
has the form (1,1,..., 1)'. Second, all non-zero eigenvalues of L are strictly 
negative. These will be denoted {— A^} with k = 1, . . . N — 1. 

One may work with the system (|1.6|) by taking a discrete Fourier trans- 
form, however we will use the eigendecomposition (jl.7p to define the so- 
called normal modes: z := Q _1 x. We readily see that these modes satisfy 
a non-interacting system of SDEs, 

dz(t) = Az(t)dt + adB(t) (1.8) 

where B = Q X W. Since Q is a orthogonal, the rows of Q 1 form an 
orthonormal family of vectors and it follows that B is a vector of independent 
standard Brownian motions. 

We observe that the mode zq, associated to the eigenvalue 0, is simply 
a standard Brownian motion. Recalling that the form of the eigenvector 
associated with the eigenvalue is (1, 1,1,..., 1)', it immediately follows that 
the "center of mass" of the bead-spring network x(t) := jjJ2n=i x n(t) = 
jfZo(t) is also a standard Brownian motion with diffusion coefficient a/y/N . 

However, we will see that individual particles are sub-diffusive processes 
with an exponent that depends on the details of the network. We trans- 
form back into real coordinates by multiplying (jl.8p on the left by Q. This 
recovers the SOU process form of each particle: 

AT-l 

Xn{t) = -j=B Q (t) + ^ Qk+l,n+lZk(t) (1-9) 
V k=l 

where the coefficients {qkn}k,n<N are the entries of Q and the {z k }k<N~i 
are defined by the system of SDEs 

dz k (t) = -\ k z k {t)dt + adB k (t). (1.10) 
1.3 Touchstone example: the Rouse chain 

We may now discuss the Rouse chain model in the context of SOU processes. 
The graph Qr associated with this model consists of edges x n <-> x n+ \ for 
all n = 1, . . . ,N, each with spring constant k. We also include the edge 
xi <-» xn so that the particles in the system are exchangeable. This yields 
the system of SDEs 

dx n (t) = n[x n -i(t) - x n (t)] + K[x n+ i(t) - x n (t)]dt + adW n (t), (1.11) 
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which can be summarized by the vector equation 



x(t) = «Lx(i) + adW(t) 
where L is the tridiagonal matrix 



V 
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The eigenvalues of the matrix L are given by 



-4sin z 



kir 

~W 



for k = 0, . . . , N — 1 (see Section 13. 4p . The fact that the eigenvalues are 
selected from what we shall call a spectral shape function 

(f(x) = 4sin 2 (7rx) 

is an essential feature to all the models studied in this paper (see Assumption 
□ in Section E]). 

The well known observation [111 \12\ [T3] that distinguished particles in 
Rouse chains are sub-diffusive with exponent v = \ follows from noting that 
for small values of x, the shape function tp(x) behaves essentially like x 2 , a 
notion that is generalized by Assumption [2l Then one can apply Laplace's 
method to the MSD to determine the asymptotic behavior, see proof of 
Theorem [2TJ 



2 Sums of Ornstein-Uhlenbeck processes 

Define the family of processes 



N-l 



XN 



(t) = Co,Ar-Bo,Af + Ck,NZk,N{t) 



(2.12) 



k=l 



where for each (k, N) the Ornstein-Uhlenbeck processes {zk,N}k=i satisfy 
theSDEs 

dzk,N(t) = -h,NZk,N(t)dt + dB kjN (t). (2-13) 
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Dependence on N will be suppressed wherever there is no chance of ambi- 
guity. 

Our focus will be on systems, such as the Rouse model, for which the 
diffusive spectrum can be analyzed asymptotically in N. The eigenvalues of 
the Laplacian matrix associated to the Rouse connection graph converged 
to the continuous shape function 4sin 2 (7rx) and we generalize the notion as 
follows. 

Assumption 1 (Diffusive spectrum shape function). There exists a non- 
negative continuous function ip G Li([0, 1]), that is strictly postive for all 
x £ (0, 1) with ip(0) = such that 



lim sup 
N-+<x> ke{i,...,N-i} 







Continuity of ip along with the specification of the value 97(0) = assures 
that the longest relaxation time tjv, which is the inverse of the smallest 
spectral value, tends to infinity with N. As is mentioned in the Discussion 
section at the end of this paper, there are natural generalizations such as 
subdiffusion in a quadratic potential where this will not be the case. 

We will see that the behavior of the shape function near zero determines 
the most important qualitative dynamics of SOU processes. 

Assumption 2 (Spectral parameter p). The shape function ip has a Frobe- 
nius expansion fWjj 



ip{x) ~ x p a-nX™ 1 for small x 

n=0 

for some p > in the sense that for each fixed N, 

lim x~ N (<p(t) ~x p J2 a n^A = 0. 

X ^ V n=0 / 

We will see (Proposition 12. 2|) that the effect of perturbations to the 
coefficient family is subdominant to the shape of the diffusive spectrum. 
In applications of interest this happens due to averaging of the coefficients, 
which we may characterize in terms of weak convergence of measures. Define 
for x £ [0, 1] the sequence of coefficient measures, 

N-l , k , 

fi N (dx) ■=^2slx- — \4 jN (2.14) 
k=o ^ ' 

where S(x) is the Dirac (5-distribution. 



S 



Assumption 3 (Convergence of coefficients). There exists a nonnegative 
finite Radon measure fi on the interval [0, 1] such that — > weakly. 

At this point, we make a note about initial conditions. The anomalous 
behavior in the limiting process is actually the infinite extension of the 
transient dynamics of the finite N processes. In order for the sequence of 
processes {x7v(i)} to be tight, it must be true that the sequence of initial 
conditions {x/v(0)} must also be tight. It is natural to choose Xjv(0) = for 
all N, but we will further simplify by choosing vanishing initial condition 
for each of the OU processes: 2fc jv(0) = 0. We note that in most relevant 
cases it is not appropriate to simply choose each z^^n from its respective 
stationary distribution. With such a choice, as N — > oo the sum of samples 
from stationary distributions will not converge. 

2.1 Asymptotic behavior of SOU processes 

We seek to relate the structure of the shape function near zero to the asymp- 
totic anomalous diffusive exponent of an SOU process. As mentioned in the 
Introduction, for any fixed, finite N, it is expected that the MSD profile will 
have the form over the three timescale regimes. Before stating a rig- 

orous description of the dynamics in Theorem l2.lt we include some intuitive 
discussion. 

The short-timescale diffusive regime has both a mathematical and a 
physical interpretation. The mathematical intuition is that Ornstein- 
Uhlenbeck processes are locally like Brownian motions. Therefore for 
t < t\ := A" 1 ^,, the process xjv(i) is essentially a finite sum of Brown- 
ian motions. Physically, in the context of distinguished particle dynamics, 
the initial diffusive regime results from the fact that for a short period the 
beads are able to diffuse independent of the constraints from the network. 

To explain the diffusive behavior on the largest timescale, we first note 
that the sum ^2^=1 c k z k converges to a stationary distribution which is 
normal with mean zero and variance X^fc^i 1 c jt/(2Afc). The timescale of the 
approach to stationarity is dictated by the longest relaxation time tn = 
\nin- F° r t > T N, the process Xjv(i) is a Brownian motion plus a stationary 
correction and so the MSD must satisfy limt_ >00 E[arjy(t)] jt = c\. 

We cannot analyze the intermediate regime exactly, but suppose that 
Assumption [2] holds for p > 0. Then \ m i n will be roughly (k/N) p , which 
means that the longest relaxation time tn = N p approaches infinity as ./V 
increases. As a result, the anomalous stage of the diffusion is increasingly 
prolonged, and is infinite in extent in the large limit. In Theorem 12. II we 



9 



show that this limit exists and in the course of the proof demonstrate that 
the MSD of the xjv(i) processes converge uniformly on compact sets. By 
performing an asymptotic analysis on the limiting process we discover the 
anomalous exponent and from the uniform convergence of the MSDs, we see 
that this is indeed the anomalous exponent seen in the intermediate phase 
of the finite N processes. 

We are ready to state the main theorem for SOU processes. For the 
large-i asymptotic statements, we say that f(t) ~ t v if lim^oo f(t)/t u = C 
for some nonnegative constant C. 

Theorem 2.1. Let the sets of processes {xiy(t)} and {zk Nif)} be defined by 
(|2.12p and (|2.13|) . respectively. Take z^ jv(0) = for all k and N . Suppose 
that the diffusive spectrum {Xkn} converges to a shape function <p in the 
sense of Assumption^ Furthermore, suppose the coefficients {c^ at} satisfy 
Assumption^ with limiting weight measure p. 

Then the family {x7v(i)} converges in distribution as N —* oo to a mean 
zero Gaussian process x(t) defined by its auto- correlation function 

E[x(t)x(s)] = f 1 6 ( 1 - e" 2 ^)^)) ^dx). (2.15) 

Jo 2v?(x) V J 

If the shape function (p furthermore satisfies Assumption with spectral 
parameter p > and the limiting weight measure fi is Lebesgue measure, 
then asymptotically, the limiting MSD function o~(t) := E[x 2 (£)] satisfies 



a(t) ~ t, t near zero. 



and 



t X ~p p > 1 

<r{t)~{ int p = i ilarge. (2.16) 

1 < p < 1 

Remark 2.1. It is important to note that the limits with respect to N and 
t implicit in (|2.16p are not interchangeable. For any finite N, K[x 2 N (t)] ~ t 
for large t because the Brownian term eventually dominates the dynamics. 
In the absence of the Brownian term (co,jv = 0), the process £/v(i) is positive 
recurrent. 

Proof. Convergence in distribution follows from establishing two standard 
facts [20j: convergence of the finite-dimensional distributions, and tightness 
in the space C([0, T]) of the family of processess {xjy} for any T > 0. Since 
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each of the processes in this sequence is Gaussian, convergence of finite- 
dimensional distributions follows from pointwise convergence of the ACFs, 
<rjv(i, s) := ^[xN(t)x]y(s)]. In order to establish tightness we will use the 
Kolmogorov criterion (|2.18p . Subsequently, the asymptotic analysis reduces 
to an application of Laplace's method [TO] . 

Convergence of the finite- dimensional distributions: We compute the ACF 
for xn. The explicit solution of the respective OU processes is given by 

z k (t) = e-^z fc (0) + f e- x ^'^dW k (t') 
Jo 

where we have suppressed the dependence of the coefficients {c^} and diffu- 
sive spectrum {A^} on N. Since the modes are assumed to be independent 
with vanishing initial conditions, we see that for s,t > 0. 

E[z k (t) Zj (s)} = 8 kj ^-e- x ^- s \ (l - e - 2A *(* As )) . 

where 5 k j is the Kronecker delta-function. Observing that cross-terms dis- 
appear and including the leading term K[Bo(t)Bo(s)] = t A s, yields 

JV-l 2 

a N (t, s) = c 2 (t As)+J2 ^e- A *l*- s l (l - e - 2X ^ tAs) ) . (2.17) 
k=l Xk 

In light of the assumption that </?(0) = 0, the above can be rewritten in 
terms of the coefficient measures defined in Assumption O 

Note that the integrand is continuous for all x E (0, 1] and can be extended 
analytically to include x = for each choice of t and s. The integrand is 
bounded above by t A s and therefore the weak convergence of the measures 
implies the limiting expression (|2.15j) . 

Tightness: As mentioned, tightness of the family of processes {xjy} in 
C([0,T]) is implied by the Kolmogorov criterion: given T > 0, there ex- 
ists anJVo 6N and strictly positive constants a, j3 and C such that 

sup E[\x N (t) - x N (s)\ a ] < C\t - s\ 1+lS (2.18) 

N>N 

for all s, t e [0, T] 
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From (|2.17p . we compute 

E[(z k (t) - z k (s))(z 3 (t) - Zj (s))] = S jk ^-(2 - e ~ 2X ^ s )(l - e' 2 ^^), 

l\ k 

while E[(i?o(i) — Bq(s)) 2 ~\ = \t — s\. Cross-terms vanish and we find 

N-i 2 

E[(x N (t) - x N (s)) 2 ] = c 2 \t - s\ + ^-(2 " e~ 2A ^))(l - e^-l) 



k=l 

N-1 2 

< cgjt _ s | + £fc(i_ e -2A*i*-.i) 
fi A * 

In the last line we applied the naive estimate (1 — e xt ) < Xt to each term 
of the sum. The sum appearing in the last line is exactly ^fc^o 1 ~ 
■^)lr 0i i](x)c| = /ijv([0, 1])- By Assumption [3l /an — ► \i weakly and by an 
equivalent statement we have 

limsup/Mr([0,l]) < /i([0,l]). 

JV^oo 

As such, there exists an No such that for all N > Nq, /xjv([0, 1]) < 1 + 
fi([0, 1]). Therefore, for all N > N , 

E[(x N (t) - x N (s)) 2 ] < (1 + /x([0, l]))\t - s\. 

Finally, noting that xjsr(t) — xn(s) is Gaussian, we see that 

sup E[{x N (t) - x N {s)) 4 } < 3(1 + fi([0, l})) 2 \t - s\ 2 , (2.19) 

N>N 

which confirms (12.180 . 

Asymptotic analysis: We now consider the large-t asymptotic behavior of 
the limiting MSD function a(t) := E[x 2 (t)] in the presence of Assumption 
[2] with shape parameter p. First we observe that for a given constant A 

A- 1 (l-e- A *)= f\- Xs ds. 
Jo 

Applying this identity to the integrand in (12.151) and subsequently using 
Fubini's Theorem to interchange the integrals yields 

a(t) = f C e- 2v{x)s p{alx)ds. 
Jo Jo 
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Recalling our assumption that fx is simply Lebesgue measure, we write 

a(t) = [ <Z>(s)ds 
Jo 

where <3>(s) is the Laplace integral &(s) := Jq e~ 2ip ^ s dx. 

Because (p(x) is continuous (by Assumption [I]) and therefore bounded, 
it follows that <&(s) is also continuous. By the Fundamental Theorem of 
Calculus, 

1 f 4 

lim- / <S>(s)ds = $(0) = 1. 
t ~^° t Jo 

This limit is finite and nonzero, which directly implies that near zero a(t) ~ 
t. 

In order to characterize large-t behavior, we note that the minimum 
value of ip{x) is assumed to be at x = and therefore the only significant 
contribution to the large s asymptotics will be in a small neighborhood near 
zero. Following |19j . for example, we have 

/ e~^ s dx~ / e- x " s dx = s~~p-Y(-) 
Jo Jo P P 

where in the last equality we applied the substitution y = x p s and T is 
the Gamma function T(z) := J °° y z ~ 1 e~ z dz. Integrating this asymptotic 
expression while minding the various ranges of values of p yields (|2.16)) . □ 

By generalizing (|2.19p to higher and higher moments, one can show that 
the limiting process x(t) is a-H61der continuous for any a S (0,1/2). This 
reinforces the notion that the limiting process is locally like Brownian mo- 
tion, but asymptotically like fractional Brownian motion. 

2.1.1 Robustness of the anomalous exponent with respect to per- 
turbing the coefficients 

As noted in [10], if one were interested in conducting statistical inference 
on a set of data, trying to fit to an SOU process, it would be a prohibitive 
task to fit the tens of thousands of coefficients. And Assumption [3] and the 
calculation for the ACF of distinguished particles processes (see proof of 
Theorem 13. 1|) can leave the false impression that the ability to calculate the 
anomalous exponent v is restricted to special, delicately balanced coefficient 
families. In fact, the result is more robust than this. We demonstrate 
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this phenomenon with the following proposition which presumes random 
coefficients. 

In what follows there are two senses of averaging. There is the probability 
space we use to select the coefficient family and for a fixed coefficient family, 
there is the probability space associated with the family of Brownian motions 
that drive the dynamics. We will use E c and Var c to denote the average and 
variance with respect to the coefficient probability space, and will continue 
to use E and Var to denote averages and variances with respect to instances 
of the Brownian motions. 



Proposition 2.2. Suppose the triangular array of coefficients {ck,N}k<N 
are independent random variables of the form c^^n = c/yN where c is a 
random variable with a finite fourth moment. Then the associated coefficient 
measures {^n}n = \ converge weakly to E[c 2 ] times Lebesgue measure almost 
surely. 

Furthermore, the conclusions of Theorem \2.1\ apply to the limiting pro- 
cess, x(t). 

Proof. The result follows from the Strong Law of Large Numbers. 

Let / : [0, 1] — ► R be a continuous function on (0, 1) and denote m n := 
E c [c n ], for n = {1,2,3,4}. We also recall the definition of the coefficient 
measures [In from Assumption (3) Then the random variables 



r 1 ( h \ 

/ f(x) m (dx) = J2f 
J ° k=l ^ ' 



c k,N 



have mean E c [Ijv] = m,2 Ylk=o f{~w) W anc ^ variance 



Var c (I 



N) 




For the second equality, we note that cross-terms of the sum vanish due 
to independence of the coefficients. It remains to recognize the Riemann 
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approximation ^2^ = q f 2 (k/N)l/N = f f 2 (x)dx + 0(1/N), which implies 
that ^ 

Var c (/ W ) = (m 4 - m|)i (J f 2 (x)dx + 0(1/N)\ , 

which tends to as N — > oo. Therefore the sequence of random variables J/v 
converges almost surely and we conclude that, in the language of Assump- 
tion El the coefficient measures fiN converge to ni2 times Lebesgue measure 
almost surely. □ 

3 Anomalous diffusion of distinguished particles 
in bead-spring networks 

3.1 General framework and main theorem 

When considering the large N scaling limit of bead-spring systems, there 
are two distinct constructions. In the polymer physics community [11] [13] . 
it is typical to simply add a new bead to the bead-spring loop while keeping 
the spring constants the same. In mathematical developments, (see [21] 
for example), it is typical to also increase the spring forces while rescaling 
magnitude of the noise in order to develop a continuum limit of the full 
bead-spring system. This is the so-called "random string" model. 

While the distinguished particle limit exists in both cases, there is a 
marked qualitative difference in the behavior of the limiting process. In the 
physics development, the limiting process is of the type described in the pre- 
ceding section: locally like a Brownian motion, but globally sub-diffusive. In 
contrast, the limiting distinguished particle process in the "random string" 
development is sub-diffusive on the shortest time scales and approaches a 
stationary distribution. As is further discussed in Section 13.61 the effect of 
"bringing the anomalous diffusion to the local scale" is that the resulting 
process is locally rougher than Brownian motion, having infinite quadratic 
variation, but finite quartic variation [22J. Our focus is not on this develop- 
ment however, because we are interesting in processes which exhibit anoma- 
lous diffusion over arbitrarily large time scales. One may be concerned that 
the full chain does not have a limit in the polymer physics construction, but 
such a limit is not our goal. Rather, we proceed with the knowledge that 
there are a finite but large number of beads in the relevant physical systems 
and although there is no convergence of the full chain structure, there is 
convergence of the effect of the chain on the distinguished particles. 
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Convergence of diffusive spectrum is sufficient for convergence in distri- 
bution of a family of SOU processes. We now argue the same principle holds 
for distinguished particle processes. The central insight is that the diffusive 
spectrum in this setting is given by the eigenvalues of the Laplacian matrix 
associated with the weighted connection graph Q. A well-established charac- 
terization of stability of a family of graphs is convergence of the eigenvalues 
to some shape function, as seen for the Rouse chain model in Section 11.31 

The only technical detail is that the coefficient structure of the dis- 
tinguished particle processes defined in Section 11.21 by Equations (II. 9p and 
(jl.lOP may not precisely satisfy Assumption [3] on the coefficient family. How- 
ever, after computing the ACF of the distinguished particle process, we see 
that it is equivalent in law to an effective SOU process to which Theorem 
12.11 does apply. 

Theorem 3.1. Let {Gn} N&§ be a sequence of graphs, each having N vertices 
respectively, with edge weight sets {£jv}jvgN such that the triangular family 
{^k,N}k<N of eigenvalues of the associated Laplacian matrices Ljy satisfy 
Assumption^ Furthermore, suppose that the graphs are constructed in such 
a way that the individual particle processes are exchangeable. 

Then the conclusions of Theorem \2.1\ hold for the family of processes 
{xo,w(i)}wGN defined by (|2.12p and (|2.13p . 

Proof. Following the notation of Section [L2l and suppressing dependence on 
N, the path of the n-th particle in the system is given by 



where the family of OU-processes {zk,N}k=i are defined by equation (|2. 13|) . 
The coefficients {q n k}k=i are the n-th row of entries of the matrix Q, which 
we recall is the orthogonal matrix whose columns are the normalized eigen- 
vectors of the Laplacian matrix L^r. 

Because the particles are exchangeable, 



where x(t) = (x\(t), . . . ,xpf(t))' denotes the full vector of all N particles in 
the system. Its dynamics are defined by the vector SDE (jl.6p and the exact 
solution is given by Duhamel's formula 




JV-l 



fc=0 



E[x n {t)x n (s)} = — E[x(t) • x(s)] 



(3.20) 
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Again, we recall the assumption that x(0) = 0. 

In order to calculate the autocorrelation E[x(i) • x(s)], we observe: 

k=1 \Jo ) k \Jo ) k 

N / N t \ / N s 



N N N 



= EEE 4 «/ 

fe=l i=l j = l 17 u J 

1/11 fc=i j=i J 
./o 

In the last line we have used the Probenius norm: \\A\\f := J2k=i Y^f=i a kj = 
X^fc ^1 wnere {^k} is the set of eigenvalues of A. 

To complete the calculation above, we note that e L * is similar to e At , 
where A is the diagonal matrix from (|1 .7j) whose entries are the eigenvalues 
of L. Therefore the eigenvalues of e L * are exactly the entries e At , namely 
the set {e _Afet }^T 1 . Imposing the assumption that x(0) = 0, this implies 

ftAs N 

E[x(t) • x(s)] = a 2 / Y^e- x ^ t+s - 2r) dr. (3.21) 
Jo k=0 

Rearranging terms, we see that for each n, the distinguished particle x n 
is a mean-zero Gaussian process with ACF of the form found in Equation 
12.171 with the coefficients identically set to c^n = a/vN. In this way, we 
see that the distinguished particle processes indexed by N are equivalent in 
law to a family effective EOU analogues {x]y(t)} defined by 



MO = ^ 

where 

dz k ,N(t) = -X k ,NZk,N{t)dt + dB kN (t). 
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The coefficient measures /xjv, defined by (|2.14p . converge weakly to Lebesgue 
measure and the asymptotic conclusions of Theorem l2.1l applv to x n directly 
and therefore to the distinguished particle process by corollary. □ 

3.1.1 Examples 

We see from the argument above that, after collecting terms appropriately, 
the coefficient family {cj^jv} of the effective distinguished particle process 
XN(t) is given by the respective multiplicities of the eigenvalue family {A& jv}. 
Whereas the coefficient measures in the Rouse model and its generalizations 
(Section 13. 2p will all converge to Lebesgue measure, we take a moment dis- 
cuss two examples where this is not the case. 

Let /Cat denote a complete graph on N vertices. Then there are only 
two eigenvalues: 0, which has multiplicity 1, and N/(N — 1) which has 
multiplicity N — 1 [23J . The distinguished particle process associated to 
each /Cat is equivalent to an effective SOU process 

1 y/N 

where 

N 

d~z N {t) = -^—z N {t)+dB N {t). 

The coefficient measures converge to the Dirac-5 distribution centered at 
x = 1. Recalling x(0) = 0, we see that the limiting MSD is given by 
E[x 2 (i)] = 1 — e _i , i.e. a one-dimensional OU process. 

The same asymptotic behavior is observed from a system with a non- 
trivial coefficient family structure. An iV-hypercube on 2 N vertices has 
eigenvalues of the form with respective multiplicities (^) [23] . By the 
de Moivre-Laplace Theorem we have the following large- ./V characterization 
of the coefficients, 

Rewriting the associated coefficients reveals a sequence of approximate 
Dirac-(5 functions centered at x = |, 

n N i 

m (dx)*J-J25 N (x)- 
k=l 

where S N (x) = yjNe' 2N ^~^) 2 . 
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Since the eigenvalue shape function is <p{x) = 2x, the limiting MSD 
satisfies 



E[x 2 (t)] = f f e- 2xs 5(x - \)dxds = 1 
Jo Jo 2 



'0 Jo 

3.1.2 Higher-dimensional diffusions 

The assumption that the particles are diffusing in one-dimensional space is 
not essential. Much like for standard Brownian motion, the dimension of 
the diffusion only has an effect on the ACF by a multiplicative constant. 
The exponent of the diffusion is unchanged and this principle holds for 
distinguished particle processes. 

Consider the family x(t) = {xi(£),X2(i), . . . ,Xjv(i)}, where for i = 
1,...,N, the particles have the form x, = [x},xf, . . . , xf) where d is the 
dimension of the space in which the particles are moving. The system-wide 
configuration potential is given by 



*(x) = - ^ K nml x 



2 



inn ■•"•-I: ■ 



In this case the interactions decouple in the various components of the dif- 
fusion. The SDE for the a-component of the n-th bead is given by 

d 

m^n (3=1 

where 5 a g is the Kronecker 5-function. We see that each component conducts 
its own diffusion independent of the other components and conclude that, 
assuming vanishing initial conditions as usual, 

2 j ft As N 

E[x n (t) • X„( S )] = ^ / V e -2A fe (t+ S -2r) dr 
iV Jo k=l 

which is simply d times the autocorrelation for one-dimensional distin- 
guished particles, as in (|3.21[) . In contrast to this observation, the dimension 
of the underlying connection graph does have significant effect, which we in- 
vestigate in Section 13.31 



3.2 The Rouse polymer model 

We return to the touchstone example from Section 11.31 t ne Rouse chain 
model. We recall that the connection graph Qr consists of edges x n «-> x n+ \ 
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for all n = {0, . . . N— 1} as well as the edge xo <-> x n . This yields the system 
of SDEs 

dx n (t) = K[x n -i(t) - x n {t)} + K,[x n+ i(t) - x n {t)} + adW n (t) 

and the diffusive spectrum is given by A^jv := 4sin 2 (^). In terms of 
Assumption [Q the spectral shape function is 

<f(x) = 4sin 2 (7rx) 

with x E [0, 1). Using the Taylor expansion for sin(x), we see that this shape 
function essentially satisfies Assumption [2] with shape parameter p = 2. 
(The only sense in which it does not satisfy the parameter assumption is 
that it is not strictly increasing on the full interval [0,1). However, this is 
easily overcome by restricting to the interval [0,1/2], and multiplying the 
MSD by two.) 

By Theorem 13. 1[ the family of distinguished processes sco,iv(*) are tight 
and the limiting process x(t) has the MSD 

E[x 2 (t)] = 2a 2 f [ 2 e- 4Ksin2( - nx >dxds (3.22) 
Jo Jo 

It follows from the conclusions of Theorem 12.11 that the anomalous exponent 
is v = ^. For an explicit development of the above Laplace integral, includ- 
ing the order of the correction terms, Eq. 13.221 is happens to be a worked 
example in [19], Chapter 6. 

3.3 Higher-dimensional Rouse analogues 

While we showed in Section [3.1.21 that large-i anomalous exponents do not 
depend on the dimension of the space in which the particles reside, we now 
observe that behavior does change if the bead-spring network has a higher 
dimensional connection graph. We employ a standard technique from graph 
theory of constructing complex graphs from simple ones [24] . 

The Cartesian product of two graphs Q\ and £/2 5 which have vertices 
{vi}f =l and {wj}*£ l5 respectively, consists of vertices enumerated by the set 
of pairs {(vi,Wj)}. There is an edge (vi,wj) <-> (vk,we) if and only if either 
Vi = Vk and Qi contains the edge Wj <-> W£, or Wj = W£ and Q\ contains the 
edge Vi <-> Vk. When applied to a cycle graph such as the Rouse graph Qr 
to itself, the Cartesian product yields the skeleton of a torus. 

The adjacency matrix of a Cartesian product of two graphs is given by 
the Kronecker sum of their respective adjacency matrices A± © Ai- If Q\ 
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is a graph with N vertices and Q2 has M vertices, this sum is defined by 
A% © A2 := Ax (8) Im + /at ® A2, where the Kronecker product A B is a 
block matrix whose blocks are of the form aijB. The resulting Kronecker 
sum matrix has dimension NM x NM. The only fact we will use here is that 
the set of eigenvalues of the Laplacian matrix associated to the Kronecker 
sum is given by the set [23] 

a%j = {A; + n s : i € 1 . . . ,N,j € 1, . . . M}. 

where {A,} and {fij} are the eigenvalues of the Laplace matrices for Qx and 
Q2, respectively. 

Denoting the iV-bead Rouse eigenvalues by {A^jv}, the MSD of a distin- 
guished particle in an N x N system is computed to be 

1 r* 

E[x N (t) 2 ] =jj2 J2 e ~ K{Xi ' N+Xj ' ds 

i,j 

which converges as N — > 00 to the integral 

E[x(t) 2 ] = f f 1 f 1 e- 4K( - sin2(7Tx)+sin2{wy))s dxdyds. (3.23) 
Jo Jo Jo 

The family of distinguished particles converges in distribution in the sense of 
Theorem 13.11 and it remains only to perform the Laplace integral asymptotic 
analysis on (|3.23jl . 

As before, we take the leading order behavior near the (x, y)-origin. For 
large s, we have 

1 1 

$( S ):=4 [ 2 [ 2 e -M^ 2 (7rx)+sm 2 (ny))s dxdy 

Jo Jo 

poo poo 

~ 4 / / e~* Kn ( x2+y2 >dxdy 
Jo Jo 

which after a conversion to polar coordinates yields that for s large, 3>(s) ~ 
s . Integrating this asymptotic relation implies that for large t, 

E[x(t) 2 ] ~ln(t). 

The Kronecker sum can be iterated arbitrarily many times for higher 
dimensional connectivity. Denoting the number of Kronecker sums by D, 
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we see that the associated <£>(s) satisfies 



$(s) ~ 2 D f exp f-4K7r 2 sV x?l dx 
+ \ £1 / 



/ e -4KTT 2 sr 2 r D-l 



dr 



10 

D 

Cs~- 



where we have allowed the constant C = C{D) to change from line to line. 
After integrating to get the MSD, we see that for D > 3, the Rouse model 
with .D-dimensional connectivity has MSD that is bounded for all time. 



3.4 Universality of the Rouse exponent 

Returning our focus to the 1-D chain, we address an observation [TT] that 
the t? long-term behavior is universal for a class of models. Such a class 
is never precisely described in the physics literature, but we provide in this 
section one characterization. 

We construct a network by starting with a Rouse chain where nearest 
neighbor edges x n «-► x n+ \ are weighted by a single spring constant k\ > 0. 



We generalize the model by allowing edges of the form x r 



n+j 



which 



are respectively given uniform weights kj > 0. The weights are assigned 
uniformly in the spirit of preserving exchangeability of the beads. After the 
appropriate edges are added at the boundaries (e.g. the edge xq <-> xn-i 
with weight K2) the resulting Laplace matrix is a circulant matrix. For 



example, under the assumption that Kj 
Laplace matrix: 



for all j > 3, we have the 
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(3.24) 



and kq = — Y^j K j ■ By convention the set of indices for the 



where n-j = 

weights {kj} will be j G {-[N/2\ , . . . , [iV/2] }. 
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Theorem 3.2 (Universality of the Rouse exponent). Let {Gn} be a family 
of graphs whose associated Laplacian matrices are circulant with weights 
{/%}jez which satisfy K_j = Kj for all j > 1 and kq = Ylj^z K j- We assume 
that there exists an integer K > 0, such that Kj = for all j > K . 

Then the family of distinguished particle processes xjy converges in dis- 
tribution to a mean zero Gaussian process x with ACF given by ()2. 15j) where 
the shape function cp is defined by 

oo 

y,( x ) : = e 2nixj Kj (3.25) 

j=-oo 

Furthermore MSD the limiting process x satisfies the asymptotic relation- 
ships 

E[x 2 (t)] ~ t, t near zero; 
E[x 2 (t)] ~ t large. 

Proof. The eigenvalues and eigenvectors of circulant matrices can be com- 
puted directly [25]. For any given set of weights {Kj}, , the associated cir- 
culant matrix has a set of eigenvalue-eigenvector pairs {(A^, ^k)}k=o given 
by 

[N/2] 



e 2mkj/(N+l) K 



3 



j=-YN/2\ 



and 



^ ._ x /-i J2km/(N+l) 2kTv2i/(N+l) 2kniN / (N +1) 

Uniform convergence of the eigenvalues to the shape function 13.251 in the 
sense of Assumption [1] is clear. Tightness of the associated distinguished 
particles processes follow as in Theorem 13.11 The near zero asymptotic 
behavior of the MSD is as in Theorem 12.11 It remains only to demonstrate 
the large-t MSD behavior. 

The restrictions that K-j = Kj and kq is the sum of all the other weights 
implies that 

K K K 

(p(x) = £ e 2mxj Kj = £ 2k j (1 - cos(27rxj)) = 4£ K j sin 2 (7rxj) 

0=-K j=l 3=1 
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We immediately see that the Rouse spectrum of Section 13.21 corresponds to 
the special case K = 1, K\ = 1. Since the sum is finite, applying the Taylor 
expansion term-by-term to the series yields 




for x < 1/K. As such, the shape function satisfies Assumption [2] with 



Readers familiar with the theory of Toeplitz matrices will recognize this 
type of weak convergence for eigenvalues from Szego's Theorem. In this 
more general light, we see that this notion of eigenvalue convergence is more 
robust than the case stated here, however if one wishes to carry out the 
program for non-exchangeable bead-spring systems, the convergence of the 
eigenvectors will have to be more carefully considered. 

3.5 The inclusion of repulsive forces 

Seeing this universal nature of the Rouse scaling discourages the notion 
that distinguished particle processes will be able to address the wide range 
of behaviors seen experimentally. It is interesting to note, at least from 
a mathematical point of view, that if we are allowed to include repulsive 
potentials between the beads, a larger class of exponents becomes available. 

For example, suppose that for a generalized Rouse network with connec- 
tivity order 2, we have K2 = — Ki/4. Then the leading order term of the 
expansion for the shape function <p vanishes. One finds that in this case, 
the spectral parameter is p = 4, and the resulting large-i MSD exponent is 



for any n G N, by Fourier inverting the shape function (p{x) = sin 2n (x). 
Such a shape function would satisfy Assumption [2] with parameter p = 2n. 
To be specific about the coefficients, we let 



for j G {-2n, . . . 2n}. It follows that Kj = ( 2 ™) . The choice to take 

only even powers of sine is made to ensure the symmetry kj = 



shape parameter p = 2. The large-i anomalous exponent v 
immediately. 



2 follows 



□ 
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3.6 Rescaled distinguished particle processes and linear 
SPDE 

As a concluding note, we return to our discussion random string model which 
leads to another set of SOU processes which have qualitatively different 
behavior. 

Consider again the Rouse chain model, but now we suppose that as the 
number of beads increases, we simultaneously increase the spring strength 
by a factor of iV 2 . Without also increasing the fluctuation strength by a 
factor of V~N, the limiting structure would collapse to a single point. That 
is, we define x n to satisfy the SDE 

dx n {t) = KN 2 {x n+1 - 2x n {t) + x n _i(t)) + aVNdW n {t) (3.26) 

In [21] the author showed that there is a non-trivial limiting object under 
the above rescaling. We define a family of functions by 

u N (k/N,t) := x k>N (t) 

and take the linear interpolation for the value of UN(y,t), for all y £ 
(lb at) - Then {u^iy, t)}^ =1 forms a tight family of functions and the 
limiting object u(y,t) satisfies the stochastic heat equation 

d t u N (y, t) = Au(y, t) + W(dy, dt) 

where W(dx,dt) is a space-time white noise [26]. The second-derivative in 
space can be anticipated by seeing the double-difference spring operator in 
(|3.26p as a discrete approximation to the Laplacian A. 

The limiting distinguished particle process is u(y, •). The exact solution 
can be expressed in the Fourier inversion 

u(y,t)= / V e^ ik y/ N e-^ 2k2 ^l N2 adB k {s) 

where we recognize this is as limit of SOU processes given by (jl.2p and (|1.3|) 
with coefficients and spectrum 

f2nk\ 2 (2M\ 
\ k ,N={- W ) , c fcjiV = co S ^— j 

One can show that these monomer paths exhibit anomalous diffusion, 
E[u(y, t) 2 ] ~ i 1 / 2 for t near zero rather than for large t. In fact, the process 
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approaches a stationary distribution for large times. As mentioned earlier in 
this section, this process u(y, •) has also received attention recently because 
its sample paths are locally rougher than Brownian motion, having finite 
quartic variation [22J. 

Returning to the discussion in the previous subsection, setting k 2 = 
— Ki/4 results in the system of SDEs 

dx n (t) = K(x n+1 - 2x n (t) + x n -i(t)) + adW n (t) 

Presumably by rescaling the spring constants by N 4 while strengthening 
the noise appropriately, one obtains the stochastic beam equation in the 
limit. One suspects that the local behavior will be rougher still and have 
anomalous exponent v = 3/4. Similar results should exist for any even 
number of spatial derivatives, but we do not pursue this line of thought 
here. 

4 Acknowledgments 

The author would like to thank Greg Forest, Jonathan Mattingly, Peter 
March and Davar Khoshnevisan for their feedback and many thought- 
provoking conversations. This work was supported in part by NSF DMS 
Grant No. 0449910. 

References 

[1] Junghae Suh, Michelle Dawson, and Justin Hanes. Real-time multiple- 
particle tracking: applications to drug and gene delivery. Advanced 
Drug Delivery Reviews, 57:63-78, 2005. 

[2] Hirotoshi Matsui, Victoria E. Wagner, David B. Hill, Ute E. Schwab, 
Troy D. Rogers, Brian Button, Russell M. Taylor II, Richard Superfine, 
Michael Rubinstein, Barbara H. Iglewski, and Richard Boucher. A 
physical linkage between cystic fibrosis airway surface dehydration and 
pseudomonas aeruginosa biofilms. Proceedings of the National Academy 
of the Sciences, 103(48):18131-18136, 2006. 

[3] Jung Soo Suk, Junghae Suh, Samuel K. Lai, and Justin Hanes. Quan- 
tifying the intracellular transport of viral and nonviral gene vectors 
in primary neurons. Experimental Biology and Medicine, 232:461-469, 
2007. 



26 



[4] R. Morgado, F. A. Oliveira, G. G. Batrouni, and A. Hansen. Rela- 
tion between anomalous diffusion and normal diffusion in systems with 
memory. Physical Review Letters, 89(10), 2002. 

[5] Raz Kupferman. Fractional kinetics in kaczwanzig heat bath models. 
Journal of Statistical Physics, 114(1-2), 2004. 

[6] Mendeli H. Vainstein, Luciano C. Lapas, and Fernando Oliviera. 
Anomalous diffusion. Acta Physica Polonica B, 39(5): 1273, 2008. 

[7] I. Santamaria-Holek. Anomalous diffusion in microrheology: A com- 
parative study. 2008. 

[8] S. C Kou. Stochastic modeling in nanoscale biophysics: Subdiffusion 
within proteins. Annals of Applied Statistics, 2(2):501-535, 2008. 

[9] John Fricks, Lingxing Yao, Timothy C. Elston, and M. Gregory Forest. 
Time-domain methods for diffusive transport in soft matter. SIAM 
Journal of Applied Mathematics, 69(5):1277-1308, 2009. 

[10] Scott A. McKinley, Lingxing Yao, and M. Gregory Forest. Transient 
anomalous diffusion of tracer particles in soft matter. Journal of Rhe- 
ology, (6), 2009. 

[11] M. Doi and S. F. Edwards. The Theory of Polymer Physics. Oxford 
University Press, 1986. 

[12] Kurt Kremer and Gary S. Grest. Dynamics of entangled linear polymer 
melts: A molecular-dynamics simulation. J. Chem. Phys., 92(8):5057- 
5086, 1990. 

[13] Michael Rubinstein and Ralph H. Colby. Polymer Physics. Oxford 
University Press, 2003. 

[14] R. Kupferman, A. M. Stuart, J. R. Terry, and P. F. Tupper. Long- 
term behaviour of large mechanical systems with random initial data. 
Stochastics & Dynamics, 2(4):533 -, 2002. 

[15] R. Kupferman and A.M. Stuart. Fitting sde models to nonlinear kac- 
zwanzig heat bath models. Physica D: Nonlinear Phenomena, 199(3- 
4):279 - 316, 2004. 

[16] R. Zwanzig. N on- equilibrium Statistical Mechanics. Oxford University 
Press, 2001. 



27 



[17] Bruno H. Zimm. Dynamics of polymer molecules in dilute solution: 
viscoelasticity, flow birefringence and dielectric loss. The Journal of 
Chemical Physics, 24(2):269-278, 1956. 

[18] Hans Christian Ottinger and Yitzak Rabin. Diffusion equation versus 
coupled langevin equations approach to hydrodynamics of dilute poly- 
mer solutions. Journal of Rheology, 33:725-743, 1989. 

[19] C. M. Bender and S. A. Orszag. Advanced Mathematical Methods for 
Scientists and Engineers I: Asymptotic Methods and Perturbation The- 
ory. McGraw-Hill, 1978. 

[20] Daniel Revuz and Marc Yor. Continuous Martingales and Brownian 
Motion. Springer- Verlag Berlin Heidlberg, 1991. 

[21] Tadhisa Funaki. Random motions of strings and related stochastic 
evolution equations. Nagoya Math Journal, 89:129-183, 1983. 

[22] Jason Swanson. Variations of the solution to a stochastic heat equation. 
Annals of Probability, 35(6):2122-2159, 2007. 

[23] Fan Chung. Spectral Graph Theory. American Mathematical Society, 
1997. 

[24] Sridhar Mahadevan. Representation Discovery using Harmonic Analy- 
sis. Morgan and Claypool Publishers, 2008. 

[25] Robert M. Gray. Toeplitz and Circulant Matrices: A Review. Founda- 
tions and Trends in Communications and Information, now Publishers 
Inc., 2006. 

[26] John B. Walsh. An introduction to stochastic partial differential equa- 
tions. In Ecole d'ete de probabilites de Saint-Flour, XIV — 1984, volume 
1180 of Lecture Notes in Math., pages 265-439. Springer, Berlin, 1986. 



28 



